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It has been suggested that magnetic fields play a dynamically-important role in core-collapse 
explosions of massive stars. In particular, they may be important in the coUapsar scenario for 
gamma-ray bursts (GRB), where the central engine is a hyper-accreting black hole or a millisecond 
magnetar. The present paper is focussed on the magnetar scenario, with a specific emphasis on the 
interaction of the magnetar magnetosphere with the infalling stellar envelope. First, the "Pulsar-in- 
a-Cavity" problem is introduced as a paradigm for a magnetar inside a collapsing star. The basic 
set-up of this fundamental plasma-physics problem is described, outlining its main features, and 
simple estimates are derived for the evolution of the magnetic field. In the context of a collapsing 
star, it is proposed that, at first, the ram pressure of the infalling plasma acts to confine the 
magnetosphere, enabling a gradual build-up of the magnetic pressure. At some point, the growing 
magnetic pressure overtakes the (decreasing) ram pressure of the gas, resulting in a magnetically- 
driven explosion. The explosion should be highly anisotropic, as the hoop-stress of the toroidal field, 
confined by the surrounding stellar matter, coUimates the magnetically-dominated outflow into two 
beamed magnetic-tower jets. This creates a clean narrow channel for the escape of energy from 
the central engine through the star, as required for GRBs. In addition, the delayed onset of the 
coUimated-explosion phase can explain the production of large quantities of Nickel-56, as suggested 
by the GRB-Supernova connection. Finally, the prospects for numerical simulations of this scenario 
are discussed. 
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INTRODUCTION 



Long-duration gamma-ray bursts (GRBs) are believed 
to be strongly asymmetric explosions of massive stars. 
Some explosion asymmetry has also been inferred in core- 
collapse supernova (SN) explosions. The asymmetry is 
expected if the progenitor star is rapidly rotating, with 
the rotation axis providing a preferred direction for a jet- 
like outflow. 

The leading theoretical model for long-duration GRBs 
based on rapid stellar rotation the collapsar model[ll, i, 
0]. In this model the core of a massive star collapses 
to form a black hole. If the star is rapidly rotating, 
the overlying stellar material can form an accretion disk 
around the black hole. Because of differential rotation 
in the disk, the magneto-rotational instability (MRI) is 
expected to develop {e.g., Ref. Q), providing angular 
momentum transport and dissipation of gravitational en- 
ergy. At the temperatures and densities (T ^ 4 x 10^° 
K, p ~ 10^" g cm"'^) present in collapsar disks, neutrino 
emission cools the gas and the collapsing outer stellar 
core accretes at rates of ~ 0.1 Mq s~^ for times ^ 10 sec, 
long enough for the jet to break out of the star and to 
power long GRBs. This accretion disk — black hole sys- 
tem thus acts as the central engine for the GRB. The 
energy source in this case is mainly accretion. Since 



most of this energy is released via neutrinos, most models 
have relied on neutrino luminosity to power the outflow. 
However, it has also been suggested that a significant 
amount of energy may be extracted magnetically and 
that strongly magnetized jets may play an active role 
in GRBexplosions [1 i, i, 0, i, i, M lH H E, Q M 



3 [13, [11) 

In order for magnetic mechanisms to be viable, a fairly 
strong fields (of order 10^^ G) are required. Such fields, 
however, are expected to be generated by the MRI-driven 
turbulent dynamo in the disk during core-collapse [^. 
Similar processes, resulting in similar field strengths, are 
also believed to be ta king place in core-collapse super- 
novae explosions [l9l |20| For example, Akiyama & 
Wheeler have argued |l9l| that the field may reach the 
level set by the equipartition with the MRI-driven tur- 
bulence, as strong as 10^^"^^ G. However, as they point 
out, this field is mostly toroidal. The large-scale poloidal 
magnetic field that is needed may require a large-scale he- 
lical dynamo [21] and will probably be somewhat smaller 
than the toroidal field. Thus, we believe it is not unrea- 
sonable for the poloidal field at the disk surface to be a 
more modest Bd ^ 10^^ G (see also Ref. Q). 

In this paper we focus mostly on an alternative sce- 
nario, in which the central engine operating inside a col- 
lapsing star is not a black hole but a young millisec- 
ond magnetar — a neutron star (NS) with a large-scale 
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poloidal magnetic field of the order of 10^^ G. Such a 
strong magnetic field can be produced by a convectively 
driven turbulent dynamo 22, 2^ inside the young NS, or 
by the MRI-driven dynamo in the differentially rotating 
collapsing core[l^. An alternative possibility is that the 
progenitor core of about 10* km has a magnetic field of 
order 10^ G, similar to the fields observed in some white 
dwarfs. When such a highly-magnetized core collapses 
into a neutron star of 10 km radius, flux freezing leads 
to the amplification of the magnetic field to 10^^ G {e.g., 
Ref. [H). 

We see that the typical values of the magnetic field 
strength, the rotation rate, and the size of the central 
engine in the classical collapsar scenario with an accret- 
ing black hole and in the millisecond-magnetar case are 
similar. The overall electro-magnetic luminosities should 
therefore be also comparable. And indeed, for a typi- 
cal surface magnetic field — 10^^ G, a rotation rate 



10 sec , and a radius i?, = 10 km, the basic 



energetics and timescales make the millisecond magne- 
tar a plausible candidate for a GRB central engine [5|. 
Ultimately, the energy source for the explosion is the ro- 
tational energy of the neutron star and the magnetic field 
acts as the agent that extracts this energy. In particu- 
lar, the rotational energy of a millisecond-period neutron 
star is of order Eiot — 5 • 10^^ erg, more than enough to 
drive a long-duration GRB. The energy-extraction time 
scale, estimated using the usual pulsar luminosity for- 
mula L 



magn 



is then of order 100 sec. 
Thus, the millisecond-magnetar central engine is essen- 
tially similar (with the magnetic field scaled up by three 
orders of magnitude and the timescale scaled down by 
six) to the Ostriker & Gunn model 2j| for powering su- 
pernovae by the spin-down magnetic power of a rapidly 
rotating pulsar. 

The idea of a millisecond magnetar as a central engine 
for gamma-ray bursts has been first proposed indepen- 
dently by Usov 25 1 and by Duncan & Thompson [22 1. It 
has been further developed by a number of authors[5 , [H, 
HHHHim,!!!,!!!]. in addition, several magnetic 
explosion mechanisms have been exp lored in the super- 
nova context [11 ll^, [2l|, 0, [sii S [11 [li, [si, [13, S 
Our present paper is also devoted to investigating the 
millisecond-magnetar scenario for GRBs, but, impor- 
tantly, viewed within the overall context of a collapsing 
star. 

First we would like to stress that the plausible overall 
energetics and timescales are, by themselves, not suffi- 
cient for making a good GRB central engine model. The 
central engine also needs to be able to produce an en- 
ergetic outflow that is (1) ultra-relativistic; (2) highly- 
coUimated; and (3) baryon-free. Most of the previ- 
ous studies have focused mostly on the energetics and 
timescales, but not on the mechanisms for producing 
an outflow that satisfies these requirements (see, how- 
ever, Refs. [1^ and 37 1 for a discussion of coUimation) . 



Also, most of these models, with the notable exception 
of Refs. [1^ and [i^, have considered a magnetar in iso- 
lation, ignoring the effect of the surrounding stellar gas 
on shaping the outflow. 

One of the major points of our paper is that the in- 
falling stellar gas is still present during the explosion and 
needs to be taken into account. Thus, an important new 
element distinguishing our model from previous works is 
the consideration of the interaction between a newly-born 
magnetar and the stellar plasma in which it is initially 
embedded. Speciflcally, we argue that the pressure and 
inertia of the surrounding gas play a key role in regulat- 
ing the magnetic extraction of rotational energy from the 
magnetar. They also force the Poynting-flux-dominated 
outflow into two collimated jets, similar to Lynden-Bell's 
magnetic towers 41 , i^l • 

In Section we discuss a specific situation relevant to 
the core collapse of a massive star: a cavity inside the 
stalled bounce-shock. The radius of such shock stays 
roughly constant allowing the magnetic fields in the cav- 
ity to grow. At the same time, both the ram pressure 
of the infalling gas and the neutrino energy deposition 
inside the cavity decrease with time. We therefore argue 
that at some point, a fraction of a second after bounce, 
the magnetic field will start to dominate the force bal- 
ance, leading to a magnetically-driven explosion. 

In order to illustrate these ideas, we introduce 
the "Pulsar-in-a-Cavity" problem as a basic-physics 
paradigm for this scenario. We describe this idealized 
problem in detail in Section . We first give a general 
description of the problem and its several versions. In 
Section , we consider the simplest case of a force-free 
rotating magnetosphere inside a fixed rigid cavity. We 
demonstrate that differential rotation of the magnetic 
field lines is inevitably established inside the cavity, even 
though the pulsar itself is rotating uniformly; as a result, 
a strong toroidal magnetic field gradually builds up. We 
then study the long-term evolution of the field inside the 
cavity and show t hat the magnetic luminosity increases 
with time. We also show that a massive, non-force-free 
plasma strip unavoidably arises in the equatorial plane 
beyond the light cylinder. We expect this phase to last 
until either the development of the kink instability mod- 
ifies the situation or until the cavity walls yield to the 
internal magnetic stresses. In Section we discuss hoop- 
stress coUimation and argue that external confinement 
and differential rotation are two important ingredients 
for coUimating relativistic Poynting-flux dominated out- 
flows. In either case, a magnetic tower forms. 

We investigate the propagation of a magnetic tower 
through the star in Section . Section is devoted to 
a general description of Lynden-Bell's magnetic tower 
model 4l| , whereas Sections — describe our modiflca- 
tion of this model for the stellar environment. Specifi- 
cally, we suggest that the external confining pressure in- 
voked by Lynden-Bell's model is provided by the gas in- 
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side a hot cocoon behind a strong shock that the rapidly- 
growing magnetic tower drives into the unperturbed stel- 
lar envelope. 

In Section , we further explore some of the 
astrophysically-interesting aspects of our model. In Sec- 
tion , we consider the transition of the magnetic tower 
expansion to the relativistic regime and the final open- 
ing angle of the tower. In Section , we discuss a possi- 
ble small-scale substructure of the magnetic tower, rep- 
resented as a "train of plasmoids" , that may be the out- 
come of the small-scale magnetic structure at the base of 
the outflow {e.g., in the magnetized corona of the mag- 
netar or of the accretion disk); it may also develop later 
as a result of MHD instabilities and flux conversion in 
the growing magnetic tower itself. In Section , we dis- 
cuss these stability issues, especially in regard to the 
Rayleigh-Taylor and kink instabilities. In Section we 
investigate the prospects for reconnection in the magne- 
tar magnetosphere or in the magnetic tower, and argue 
that reconnection is ineffective in the dense environment 
deep inside the collapsing star. In Section , we address 
an important issue of ^^Ni production and argue that the 
two-phase nature of the explosion in our model is well- 
suited to explain a large amount of ^^Ni inferred from 
observations. In Section we discuss the implications of 
our model for pulsar kicks. Finally, in Section , we sug- 
gest some directions for future numerical simulations of 
this problem. We draw our conclusions in Section . 



MAGNETAR INSIDE A COLLAPSING STAR: AN 
OUTLINE OF THE GENERAL SCENARIO 

Numerous studies of core-collapse supernovae have 
shown that, as the core of a massive star collapses into 
a proto- neutron star (PNS), a bounce shock is launched 
back into the star but quickly stalls at about 200 km (see, 
e.g., Refs. (isl . |4^). The explosion then enters a long 
1 sec) phase (see Fig. [1]) during which the shock is 
quasi-stationary as determined by the balance between 
the ram pressure of the infalling material, which tends 
to quench the shock, and the thermal pressure of the 
post-shocked gas, supported by the continuous neutrino 
heating. Eventually, if neutrinos win, the shock engulfs 
the entire star and one gets a successful SN explosion. 
If they lose, the shock dies and the PNS collapses into 
a black hole that subsequently swallows the rest of the 
star, without a SN. 

In our model, we add a third dynamical component to 
this picture — the magnetic field. The magnetic pres- 
sure is pushing out, helping the explosion, as is the ther- 
mal pressure of the neutrino-heated gas. Our main idea 
is that these two outward forces evolve differently with 
time, and thus the explosion may be a two-stage pro- 
cess. In particular, we suggest that the magnetic pres- 
sure is not important during the first few hundreds of 



msec of the stalled-shock phase. However, during this 
time the magnetar makes several hundred revolutions. 
This results in a great amplification of the toroidal mag- 
netic fiux by the differential rotation (see Sec. ), whereas 
both the neutrino energy deposition and the accretion 
rate gradually decline (see Fig. [J). For example, assum- 
ing Ro = 3i?Lc = lOi?* = 100 km, and = lO^^ G, 
the entire cavity is filled (see Sec. ) with 3 • 10^^ G 
fields after about 100 turns (0.1 sec), corresponding to 
the magnetic pressure of 4- 10^^ erg/cm^. The ram pres- 
sure of the infalling stellar material a,t r = Rq, com- 
pressing the magnetosphere, can be estimated roughly 
as Pram ~ Mwff/47ri?2 ~ 8 • 10^7 Mo Mg/^ Ro^r^evg cm'^ 

where vg = (2GM/i?o)^/^ =i 5 • 10^ Mq^^ i?c;;7/^ cm/sec 
is the free-fall velocity at radius Rq, and Mq and Mq are 
the mass within Rq and the accretion rate at this ra- 
dius, expressed in units of Mq and Mq/sbc, respectively. 
This shows that, after a delay of a few hundreds of mil- 
liseconds, the magnetic pressure inevitably becomes an 
important driving force and may lead to a successful ex- 
plosion. 

To summarize our picture, the ram pressure of the ac- 
creting material provides a nurturing womb in which the 
baby magnetic field grows, until it is finally strong enough 
to break out. Neutrino energy deposition plays an im- 
portant role during this gestation period, as it prevents 
the magnetosphere from being completely squashed by 
the accreting gas. Finally, if the above picture is correct 
and the explosion does become magnetically-driven, then 
the hoop-stress mechanism makes it highly collimated, 
thus satisfying one of the key necessary conditions for 
GRB (see Sec. ). Note that this jet is driven by the 
magnetar-level (i.e., ~ 10^^ G) field and is thus much 
stronger and faster than the LeBlanc- Wilson jetjs^] that 
may have been launched a few seconds earlier, during the 
core-collapse process [30|. 



THE PULSAR-IN-A-CAVITY PROBLEM 

As a first step in trying to understand how a 
millisecond-magnetar central engine works in the coUap- 
sar context, let us consider the following basic plasma- 
physics problem: an axisymmetric pulsar inside a con- 
ducting cavity with a low-density plasma (see Fig. [3]). 
Specifically, we consider the cavity's radius i?o to be is 
much larger than the pulsar light-cylinder radius i?LC- 
We call this idealized problem the Pulsar-in- a- Cavity 
problem 3, S^l . It is a modification of the Goldreich & 
Julian model for an isolated pulsar's magnetosphere [i^; 
it has direct connections to the model proposed by 
Ostriker & Gunn for powering longer-term supernova 
lightcurves[2J| and it is also related to the models consid- 
ered by Kardashevji^ and by Illarionov & Sunyaev 47 1. 

Note that the behavior of Pulsar- in- a- Cavity depends 
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on the cavity properties. For definiteness, we assume 
that both the waUs and the plasma inside the cavity are 
perfect conductors. We also assume that all the field lines 
close back to the pulsar inside the cavity, keeping in mind 
that the magnetic field had been generated inside the NS 
and then emerged through its surface. 

At the same time, we are dealing with a whole fam- 
ily of problems distinguished by the assumed mechanical 
properties of the cavity. For simplicity, we shall concen- 
trate on the case of a spherical cavity with rigid walls (see 
Sec. ), which may correspond to early stages of the sys- 
tem's evolution. In the future we plan to consider more 
realistic but also more complicated cases where the shape 
and the size of the cavity are not fixed but instead are 
governed by a balance with an external pressure (similar 
to Lynden-Bell's magnetic tower model[4l[), and, finally, 
will consider a young pulsar in a fully dynamic environ- 
ment of a collapsing star (c.f., Ref. (24] ). 

Each of these versions of our pulsar-in-a-cavity prob- 
lem will provide important insights into the workings of a 
millisecond magnetar inside a collapsing star. They will 
probably require numerical simulations using relativistic 
force- free or relativistic MHD codes (see Sec. ). To set 
the stage for these numerical studies, we will, in this sec- 
tion, qualitatively discuss a plausible physical picture of 
the system's evolution. 

To investigate the interaction between the central 
magnetar and the surrounding stellar material, a full 
magnetohydrodynamic (MHD) description that includes 
plasma pressure and inertial effects will eventually be re- 
quired. This is especially relevant if there is a strong 
wind driven off the PNS by neutrinos and /or by the 
magneto-centrifugal mechanism {e.g., Refs. [33l.l37l.l38t). 
Of particular interest is the confinement of the expand- 
ing magnetosphere by the surrounding plasma and the 
dynamical response of the star to the expanding mag- 
netosphere at its center. For simplicity, however, here 
we shall limit ourselves to the relativistic force-free case 
in which electromagnetic forces dominate the dynamics 
almost everywhere inside the cavity (except at the the 
equatorial plane outside the light cylinder — see below) . 
While it is not realistic in the central region of a col- 
lapsing star, the force-free description may nonetheless 
refiect some essential physics, especially for late phases 
of the evolution. 

We would like to remark that the ideal-MHD assump- 
tion is well justified due to the very large plasma den- 
sities and temperatures. Specifically, the high plasma 
density ensures that the plasma is highly coUisional and 
hence is well described by resistive MHD; other non-ideal 
terms in generalized Ohm's law are unimportant. On the 
other hand, the resistivity is actually small in absolute 
terms {Rem ^ 1). AH this makes ideal MHD a good ap- 
proximation in the environment of a collapsing star[15|. 
Therefore, in this paper we shall ignore any non-ideal 
effects, leaving them for a future study. 



As a final point, in contrast with the isolated pulsar, 
which is usually regarded as stationary, our pulsar-in-a- 
cavity problem is intrinsically time-dependent. 

Pulsar in a Fixed Spherical Cavity 

In this section we consider the case of a spherical cav- 
ity whose radius, Rq, is kept constant (or is changing 
slowly). We are interested in the evolution of an ini- 
tially dipole-like magnetosphere after the pulsar is spun 
up suddenly. In the isolated pulsar case, the field lines 
extending beyond the light cylinder are swept backwards 
and openj^; there is no feedback of the outer region 
on the inner magnetosphere. In our case, however, the 
entire magnetosphere is contained inside the cavity and 
hence remain closed at all times. This is a very important 
difference between the two cases. 

The first thing to note is that a differential rotation 
is established in the magnetosphere. Indeed, consider a 
field line ^' that extends beyond the light cylinder (Fig.[3|) 
and compare the angular velocities at two points: point A 
where this line attaches to the pulsar and point B where 
it intersects the equator. The angular velocity at point A 
is clearly equal to that of the neutron star: D.a = f^*. 
Next, what is the angular velocity at point Bl Since the 
field line in question extends beyond the light cylinder, 
it cannot remain purely poloidal and a toroidal field has 
to develop. Due to the assumed symmetry with respect 
to the equator, however, the toroidal field has to van- 
ish at z = 0, B^{z = 0) = 0, and the plasma cannot 
slide toroidally backward along the field. Therefore, 
has to equal the angular velocity of the plasma at this 
point. Now, the toroidal field that develops in the mag- 
netosphere off the equatorial plane continuously brakes 
the star down, so that there is an outward flux of angular 
momentum and a Poynting flux of energy along the field 
line. The rotational energy of the pulsar extracted by the 
magnetic field is partly accumulated in the magnetic form 
and partly transferred to the equatorial plasma. Thus, 
the material at point B is continuously torqued up by the 
magnetic field. Then, since the confining wall prevents 
the material from moving out freely in the radial direc- 
tion, the toroidal velocity of the plasma becomes closer 
and closer to the speed of light. However, it can never ex- 
ceed the speed of light; therefore, the angular velocity at 
point B is bounded: — c/Rb = ^^Rlc/Rb- Thus, 
we see that the field line experiences differential rotation 
at a rate An = - ^b > fi*(l - Rlc/Rb); For field 
lines with Rb ^ RhC we get Afi « O*. This differen- 
tial rotation is established on the cavity's light-crossing 
time-scale, to = Rq/c^ 

This differential rotation is important because it leads 
to a continuous injection of toroidal magnetic flux (of 
opposite signs) into the upper and lower hemispheres. If 
the cavity's size is fixed, then the toroidal magnetic field 
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at any point grows roughly linearly with time. This is in 
sharp contrast with the unbounded pulsar case in whcih 
a steady state is established on the light travel time-scale. 

Let us now analyze the field evolution on long time 
scales {t ^ to = Rq/c) and at distances R ^ Rhc- 
Note that the toroidal magnetic field continuously in- 
creases, whereas the poloidal magnetic field docs not. 
The poloidal electric field, -Epoi, may become much 
larger than Bpoi but in any case cannot exceed the 
value Bpoin^Ro/c — i?poii?o/^LC- Thus, after several 
light-crossing times {t ^ t^) the magnetosphere out- 
side the light cylinder is dominated by the toroidal field, 

Next, since B^ ~ t, the relative change in i?^ over At ~ 
to becomes small at late times, t ^ to. An approxi- 
mate force-free equilibrium is then established in each 
hemisphere, described by the relativistic force-free Grad- 
Shafranov equation. In the i?0-dominated limit this 
equation reduces to //'(^') = 0, where / = RB^ is the 
enclosed poloidal current. The obvious solution of this 
equation is = Iq — const, i.e., a singular line cur- 

rent Io{t) along the rotation axis. The toroidal magnetic 
field is the vacuum field produced by this line current, 
B^{t, R, Z) = Io{t)/R; i.e., constant on cylinders. The 
main force balance in the magnetosphere is between the 
toroidal field tension and pressure. In other words, the 
poloidal current becomes spatially separated from the 
toroidal magnetic field: it fiows out of the pulsar along 
the axis (in both hemispheres), then as a surface current 
along the cavity walls, and finally returns to the pulsar 
along the non-force-free equatorial current sheet present 
due to the sharp reversal of the toroidal field across the 
equator. The current density in the bulk of the magne- 
tosphere is relatively small. This is similar to the electric 
current structure of the magnetic bubble considered by 
Lyutikov & Blandford[ii|. 

We shall express magnetic quantities characterizing the 
field in the cavity in terms of the total poloidal magnetic 
fiux extending beyond the light cylinder. This flux 
can be crudely estimated from the pure dipolc magnetic 
field, i.e.. 



line \E': 



~ 5'dipolc(^Lc) = ^* 



Ri_ 

Rlc 



(1) 



Then, the characteristic poloidal magnetic field strength 
in the cavity at distances r ~ i?o from the center and off 
the equatorial plane can be estimated as 



B, 
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Bo 



Ro' 



B, 
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(2) 



Let us now estimate the poloidal line current /g {t) and 
hence the characteristic toroidal field in the cavity. The 
poloidal current is found by following the shape of a field 



/(^',t) = Am 



d/poi 



^pol^^('pol) 



(3) 



where Zpoi is the path-length along the poloidal field. The 
main contribution comes from large distances, R ^ Rq, 
and thus, using Afl ~ 57» = c/Ri^c, we get 



Io(t) - 



Ro Rlc to 



(4) 



We see that for t ^ to the poloidal current becomes much 
stronger than that in the isolated pulsar magnetosphere 
(/ ^ *o/^Lc)- Using the estimate ([T]) for ^E'o, we obtain 



Io{t)^B, 



Ri t 



(5) 



Correspondingly, the characteristic toroidal magnetic 
field at distances of order Rq is 



B^{Ro) - — 

Uq 



Bo^*t , 



(6) 



which is similar to the estimate presented by Kardashev 
for the toroidal field of a pulsar inside an expanding su- 
pernova cavity [4^. Thus, after many light-crossing times 
across the cavity, B^{Ro) becomes much larger than the 
toroidal field of an isolated pulsar at these distances 
[^|oiatcd _ ^o/(i?oi?Lc) = BoiRo/RLc) = Boft.to « 
Bo^*t]. 

As we noted above, the magnetosphere outside the pul- 
sar light cylinder cannot be entirely force-free. Because 
the toroidal magnetic field reverses across the equator, 
the magnetic tension continuously accelerates the equa- 
torial plasma in the toroidal direction. The tension force 
performs mechanical work on the equatorial plasma and 
so a part of the rotational energy extracted from the pul- 
sar is deposited in the equatorial plane (the rest is stored 
magnetically in the bulk of the cavity) . Since the plasma 
in the equatorial plane rotates ultra-relativistically, the 
added energy leads to an increase in the relativistic 
"mass" of the plasma. Am ~ t'^. Correspondingly, this 
relativistically rotating massive equatorial sheet experi- 
ences an outward centrifugal force, Fcont- This force can- 
not be balanced by the toroidal magnetic field because 
the latter is zero at the equator and so the equatorial 
plasma moves out towards the wall. It then pushes the 
poloidal magnetic flux out and concentrates it a narrow 
equatorial strip of ever-decreasing width d{t) <C i?o near 
the wall (see Fig. |4]). Because of this, nearly all the 
poloidal flux \l/o that extends beyond the light cylinder 
crosses the equator at cylindrical radii i? ~ i?o. At the 
same time, in the magnetosphere above and below the 
equatorial plane, the poloidal field lines that emanate 
from this strip fan out to fill the cavity volume. Thus, 
the characteristic poloidal magnetic field in the cavity is 
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of the order Bq — "^o/Rq (see eqn.[2]) and is much weaker 
(by a factor of d/Ro) than in the equatorial strip. 

Let us assess the centrifugal force quantitatively. The 
total torque on the massive equatorial strip is r(t) = 
/ /(^',t)d\E' ~ /o(t)^o and the total work per unit time 
(the total Poynting flux coming to the strip) is Pstrip — 
Tc/Ro = BffjC'^Q. This power goes into accelerating the 
plasma rotation, that is, into increasing the rotational 
7-factor and hence the relativistic mass m of the plasma 
in the strip: 



m(t)(? 



Pt 



Rq *o ^LC R2/.Nn3 /7^ 
~R~~R -^^A*)^0- {') 



That is, the kinetic energy in the equatorial strip is al- 
ways small compared with the magnetic energy in the 
cavity, B'^{t)R^. The centrifugal force acting on the 
equatorial strip is 



Fccnt{t) 



m{t)c^ 

Rn 



BlRlnie 



R 



LC 



Ro 



Bimi 



R 



LC 



Ro 



(8) 

This force grows quadratically with time, just as the 
toroidal field pressure, but always remains small (by a 
factor of i?Lc/^o 1) compared with the total horizon- 
tal force exerted on the side wall by the toroidal field. 
However, since -Fccnt is concentrated in the thin equato- 
rial region, it may be important in a subsequent expan- 
sion of the cavity (c.f. non-relativistic MHD simulations 
by Matt et ai, Ref. Q). 



A detailed analysis of the internal structure of the mas- 
sive equatorial plasma strip is an interesting problem that 
should be studied but it lies beyond the scope of this pa- 
per. 

Another very important point is that the rate at which 
the magnetic field in a confined magnetosphere extracts 
rotational energy from the central rotating conductor ac- 
tually grows with time. This is because the magnetic 
torque per unit area is proportional to the toroidal field 
at the conductor's surface and the latter grows linearly 
with time. Thus, as long as the cavity does not expand 
(or expands slowly) and the rotation rate of the spinning 
pulsar stays constant, the magnetic power generated by 
the pulsar inside a cavity increases linearly with time: 



p{t) = /(t)i'or2, = nit 



Ro 



Pi 



isolated 



ct 

Ro 



(9) 



where -Pisoiatod B'lR^n'l/c^ is the spin-down power of 
an isolated, unbounded pulsar. Hence, after many light- 
crossing times, the power of a pulsar-in-a-cavity greatly 
exceeds that of a classical isolated pulsar. This is resolves 
the apparent paradox raised by Lyutikov 
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This runaway behavior can be attributed to a positive 
feedback between the energy that has been already ex- 
tracted from the pulsar, and the strength of the agent 
that extracts the energy (the toroidal magnetic field). 



Namely, most of the extracted energy is stored in the 
toroidal magnetic field, and, since the volume is finite, 
the toroidal field strength increases with time. Since 
the magnetosphere remains in a quasi-equilibrium, the 
toroidal field constantly readjusts everywhere and the 
inner magnetosphere feels the presence of the outer con- 
fining wall. In particular, the toroidal field at the pul- 
sar surface increases linearly with time, and so does the 
magnetic spin-down torque on the pulsar. This picture 
is similar to the combustion chamber of a rocket. In that 
case, the gas temperature and pressure increase as the 
chemical energy of the fuel is released in the combustion 
process. At the same time, the fuel burning rate grows 
with the ambient temperature. Therefore, a rapid and 
efficient burning demands high pressure and is hence fa- 
cilitated by a strong confining chamber. Similarly, in our 
case of a pulsar placed inside a cavity, the presence of 
strong cavity walls leads to an increased energy extrac- 
tion rate from the pulsar. 

In reality, we don't expect this power growth to last 
indefinitely. It may saturate, for example, due to the 
development of the kink instability, resulting in the con- 
version of the toroidal flux to poloidal flux and to the 
dissipation of some of the magnetic energy (see Sec. for 
more discussion). 



Hoop-stress collimation: contrast with the isolated 
pulsar 

The toroidal field generated by the differential rotation 
exerts a constantly-growing pressure on the cavity walls. 
If we now relax the fixed-wall assumption, this pressure 
will inflate the cavity. Will this inflation be isotropic or, 
say, collimated along the axis? 

Generally speaking, since the toroidal field pressure in 
the horizontal direction is partly negated by the field's 
tension, one expects the resulting expansion to be mostly 
vertical. However, the differential rotation producing B^ 
is relativistic: Ar2_Ro ~ ri*i?o 3> c, and it is well-known 
that hoop-stress collimation is not a trivial issue in the 
relativistic case. Thus, it is not immediately obvious this 
mechanism can be applied to our pulsar-in-a-cavity sce- 
nario. The quintessential example of this lack of colli- 
mation for ultra-relativistic magnetically-dominated out- 
flows is the isolated aligned pulsar wind inside the ter- 
mination shock. The basic reason for this is the decolli- 
mating force due to the poloidal electric field, E'poi- In- 
deed, in the case of an unbounded relativistic uniformly- 
rotating force-free magnetosphere in a steady state, i?poi 
and B^ are nearly equal at large distances from the 
axis [45}. Importantly, this balance happens in an un- 
coUimated, quasi-spherical poloidal magnetic field con- 
figuration; an excellent example of this is Michel's split- 
monopole solution [4^. Here is a crude argument explain- 
ing this lack of hoop-stress collimation in the relativistic- 
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rotation case. Consider an uncoUimated field configu- 
ration; the poloidal magnetic field is open outside the 
light cylinder and has a split-monopole geometry, i.e., 
drops off as . In a steady state, the poloidal electric 
field is i?poi = BpoiR/RhC, where R is the cylindrical 
radius; it hence drops off along radial rays as r~^. But 
also drops off as r~^. Moreover, at the light cylinder, 
-E-poi = ^poi ~ Bcj,. Since outside the light cylinder they 
both decrease as r^^, they remain comparable to each 
other at large distances. In fact, as Goldreich & Julian 
showed (is! . Epoi and B^ become equal asymptotically 
as r — > oo. The bottom line is that a quasi-spherical 
relativistic force-free equilibrium can be established as 
a balance between the collimating pinch force (the sum 
of the toroidal magnetic field pressure and its tension) 
and the opposing electric force. Hoop-stress coUimation 
is suppressed as a result of this balance. 

On the other hand, the case of a rotating magneto- 
sphere enclosed inside a cavity is different and hoop-stress 
collimation can in fact work. Indeed, as we showed above, 
at late times the toroidal magnetic field filling the cavity 
becomes stronger than both Bpoi and Spoi, in contrast 
to the isolated pulsar case. Furthermore, this toroidal 
field is distributed nonuniformly; namely, B^ R~^. 
Correspondingly, the magnetic pressure pushing verti- 
cally against the top and bottom walls is much higher 
than that on the side walls. Therefore, if we now al- 
low the cavity to expand under this pressure, the expan- 
sion will be mostly vertical. The situation is then sim- 
ilar to the non-relativistic magnetic tower proposed by 
Lynden-Bell 4l| . We therefore envision that long-term 
result will be the creation of a pair of oppositely-directed 
magnetic towers [l5j|. The interaction of the expanding 
towers with the surrounding stellar envelope aids in their 
confinement, similarly to jet collimation i n hy drodynam- 
ical simulations of the coUapsar model [3,lM|5ll. In the 
scenario considered here, the towers are driven not by 
a differentially-rotating disk, but by a rapidly-rotating 
magnetar. This suggests that considering the pulsar 
magnetosphere inside a cylindrical, as opposed to spher- 
ical, cavity may represent yet another interesting topic 
for future research. 

An important element in the above discussion is the 
fact that -Epoi is small compared with i?^. This is because 
B^ is generated as a result of differential rotation. This 
highlights the important role of differential rotation (as 
opposed to uniform relativistic rotation) in collimating 
relativistic force- free outflows. 



MAGNETIC TOWER INSIDE A STAR 

For simplicity, in this section we shall mostly con- 
sider the case when the central engine is an accretion 
disk around a black hole. However, we believe that the 
millisecond-magnetar case is essentially similar. 



Lynden-Bell's Original Magnetic Tower Model 

We suggest that the model most naturally suited to de- 
scribe the propagation of a Poynting-flux dominated jet 
through a star is the magnetic tower model, originally in- 
troduced by Lynden-Bell in the AGN context 4ll.[4a|. A 
magnetic tower is an axisymmetric magnetic configura- 
tion that arises when a system of nested closed flux sur- 
faces, anchored in a differentially-rotating disk, is twisted 
and, as a result, inflates, but when this inflation is con- 
trolled by a surrounding external pressure. The basic 
physical mechanism of this process can be described as 
follows (see Fig. [5]). 

Consider a thin conducting disk with some vertical 
magnetic flux frozen into it. Let us assume that ini- 
tially the magnetic field has a dipole-like topology (see 
Fig. [5^) , with the two footpoints of each field line located 
at different radii on the disk. Now let the disk rotate non- 
uniformly {e.g., a Keplerian disk). Then, each field line ^> 
is twisted at a rate Ai7(\I') equal to the difference in an- 
gular velocities of its two footpoints. Correspondingly, 
toroidal magnetic flux is generated from the poloidal flux. 
The pressure of the toroidal fleld pushes the flux sur- 
faces out, against the poloidal field tension. It is assumed 
that during the initial stages of this process the gas pres- 
sure, as well as the gravitational and inertial forces, are 
negligible in the magnetosphere, so the magnetic field 
is force-free. Then the expansion is uncoUimated, typ- 
ically at a 60° angle with respect to the rotation axis. 



as shown in Figure [lb (e.g.. Refs. 52, l53|). However, as 
was shown by Lynden-Bell[41[, if there is some, no mat- 
ter how small, external gas pressure Pext surrounding the 
expanding magnetosphere, then the sideways expansion 
ceases once B'^/Stt drops down to R^t- Subsequently, as 
again was shown by Lyndcn-Bell[4]J, the twisted mag- 
netosphere continues to expand in the vertical direction 
(Fig. O:) and eventually forms a cylindrical column that 
Lynden-Bell termed a magnetic tower (see Fig. [5jl). If 
the external pressure outside of the tower is kept con- 
stant and uniform, then the top of the tower rises at a 
constant speed. Plasma inertia never plays any role; the 
entire evolution is a sequence of force-free magnetostatic 
equilibria with a pressure balance between the external 
gas outside of the tower and the magnetic field inside. 

Note that the assumption that both ends of the field 
lines connect to the disk is not essential. A similar be- 
havior is encountered in the case of a rotating conduct- 
ing disk magnetically connected to a central star {e.g., 
Refs. 53, 54 1) or a rotating black holef55]. We be- 
lieve that the model is also applicable to the millisecond- 
magnetar central engine scenario for GRBs. 

To get a physical feeling of how the magnetic tower 
grows, it is instructive to derive some simple order- 
of-magnitude estimates and scaling relationships. The 
main input parameters that set the characteristic phys- 
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ical scales are the total poloidal magnetic flux "ifo (per 
unit toroidal angle) in the tower, the characteristic dif- 
ferential rotation rate Afi, and the external pressure Pext- 
Let us now ask how the main parameters of the tower, 
namely, its radius Rq, the typical magnetic field Bq, and 
the growth velocity Vtop, scale with the three input pa- 
rameters. 

First, the radius of the tower, i?o, and the characteris- 
tic poloidal magnetic field, Bpoi are related via 



pol 



rig 



(10) 



The radius adjusts so that the magnetic pressure inside 
the tower equals Poxt- From the force- free balance inside 
the tower we expect -Bpoi ^ ^o! hence, the total 

magnetic field strength at the outer edge of the tower is 
also of order Bq . Then, from the the condition of pressure 
balance across the tower's wall we get 



B. 



and so 



( 



(11) 



(12) 



To estimate the growth rate of the tower, note that the 
toroidal magnetic flux x is continuously generated from 
the poloidal flux by the differential rotation 



(13) 



Taking the tower to be a cylinder with radius 
i?o and height Ztop, we get i?^ ~ x/i?o^top = 
(^'o/i?o^top) ArJt = Bo Af2t(i?o/^top)- However, as 
stated earlier, the typical toroidal held in the tower 
should be of the order of B^; therefore, the height of 
the tower increases steadily as 



(14) 



In other words, the tower grows at the speed of order of 
the typical differential rotation velocity Rq A57. If the ex- 
ternal pressure does not change, the radius of the tower, 
determined by Eq. stays constant during its growth; 
therefore, after many turns (AOt » 1), ^top ^ Ro^ *-e-, 
the tower becomes slender. 

Since the first analytical solution proposed by Lynden- 
Bell[4l|, the magnetic tower concept is becoming more 
and more accepted by the astrophysical community. The 
formation and evolution of ma gne tic towers have been 
studied in numerical simulations 
and even in real laboratory experiments BsT 64|. 

An interesting question is the fiow of energy through 
a magnetic tower. As one can easily see, Foynting flux 
flows up from the disk along the inner segment of each 



field line and down to the disk along the outer seg- 
ment. Indeed, for each field line VE", the inner, faster- 
rotating footpoint (1) performs work on the magnetic 
field, Wi ~ Jpoi[^'(l)]ri(l) (per unit time and unit 
poloidal flux). The corresponding decelerating torque 
per unit flux is n ~ /poi[5'(l)]. In turn, the mag- 
netic field exerts an accelerating torque per unit flux 
T2 ^poi[^(2)] on the outer disk footpoint of the same 
field line. Correspondingly, it performs work at a rate 
W2 ^ /poi[^'(2)]ri(2). Because of force-free equilibrium 
in the tower, /poi[*(l)] = /poi[«'(2)] = /poi(^'), and so 
Ti = T2, i.e., all the angular momentum extracted mag- 
netically from point 1 is transferred to point 2. The two 
energy fiows, on the other hand, are not equal: since 
ri(l) > 1^(2), the energy extracted from point 1 along 
the inner segment of the field line is greater than the 
energy that fiows down along the outer segment and is 
deposited in the disk at point 2. The difference, propor- 
tional to /poi(vE')Ar2(\l/), is the power driving the expan- 
sion of the tower. A part of it goes into filling the growing 
volume of the tower with magnetic energy, and the rest 
goes into performing work against external gas pressure 
and driving the shock through the star. 

It is interesting to note that the total vertical Foynt- 
ing flux in the two segments only involves the differential 
rotation Af2 = $1(1) — $1(2), but is independent of the 
absolute rotation itself. This is because we are dealing 
here with a force-free equilibrium, so that /poi is con- 
stant along the entire length of a field line; in particular, 
it has the same sign on the two segments of the field 
line, and hence so does B^. The situation is drastically 
different in the relativistic-rotation case where both field- 
line segments extend beyond their respective light cylin- 
ders. In that case, one no longer has a force-free equilib- 
rium along the entire field line; in particular, equilibrium 
breaks down at the farthermost tip of the line where the 
two segments join. As a result, the signs of /poi (and 
hence of -B^) on the two segments are opposite, which 
corresponds to both segments being swept back. Con- 
sequently, the Foynting fiux is outward along both seg- 
ments. A similar situation arises in the non-force-free 
MHD case; the two field-line segments are then swept 
back by plasma inertia if they extend beyond the Alfven 
point. This again results in an outward Foynting flux 
along both segments. In both of these cases, the total 
vertical Foynting flux depends on the absolute rotation 
rates $1(1) and $1(2) themselves, as opposed to just their 
difference. 



Magnetic Tower Driving a Shock through a Star 

There are several reasons that make the magnetic 
tower an attractive model for the formation and prop- 
agation of a magnetically-dominated jet through a star 
within the coUapsar model for GRBs (and core-collapse 
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SNe). First, a configuration where aU the field lines close 
back onto the central engine (an accretion disk or a mag- 
netar) is natural for a field created by a dynamo with zero 
net flux. In addition, mixing of baryons from the stellar 
envelope via Kelvin-Helmholtz and/or Ray leigh- Taylor 
instabilities may be inhibited by the magnetic field. 

In order to apply the magnetic tower model to the 
coUapsar scenario, we first want to make some modifica- 
tions to Lynden-Bell's picture [3]. Specifically, we take 
into account the high-pressure cocoon that surrounds 
and confines the tower (see Fig. O. In our model, the 
magnetic tower grows very rapidly and acts as a piston 
driving a shock ahead of itself. The shocked gas above 
the tower has very high pressure; it squirts sideways and 
forms backflows that fill the cocoon around the tower. 
Therefore, the external pressure confining the tower is 
no longer an arbitrary parameter, as in Lynden-Bell's 
model, but is determined by the jump conditions across 
the shock surrounding the cocoon and across the contact 
discontinuity between the cocoon and the tower. The 
external unperturbed pressure of the star is actually ir- 
relevant; it should thus be excluded from our three input 
parameters. Instead, the expansion is controlled by the 
ram pressure related to the gas inertia; therefore, we re- 
place Pext by the unperturbed stellar density po in the list 
of basic dimensional parameters (along with '^q and Afi) 
that determine the physical scales in our problem. This 
change is an important difference between our model and 
Lynden-Bell's. 

The actual situation is complicated further by the two- 
dimensional character of the problem. Since the sound 
speed in the cocoon is very high, gas pressure tends to 
be equalized throughout the cocoon. This expectation is 
supported by th e hy drodynamic simulations of the col- 
lapsar model [s^, [Elf, which show a relatively weak (just 
a factor of 5 or 10) variation of the cocoon pressure along 
its length. This is very moderate compared to the cor- 
responding variations of the unperturbed stellar density 
and pressure, which both vary by many orders of mag- 
nitude along the vertical extent of the cocoon. Thus, 
the gas pressure is very high everywhere in the cocoon 
and so the cocoon also drives a sideways shock into the 
star. The boundary between the tower and the cocoon is 
a contact discontinuity, whereas the boundary between 
the cocoon and the rest of the star is a two-dimensional 
strong shock of some complicated shape. 



Simple Estimates 

Let us now show how the basic parameters of the grow- 
ing magnetic tower scale with \l/o, Afi, and po- We shall 
ignore any non-uniformity of the gas pressure in the co- 
coon. Also, for simplicity we shall use one-dimensional 
shock jump conditions. Since the pressure of the unper- 
turbed stellar gas upstream of the shock is neglected, the 



shock is strong. Assuming an adiabatic index of 5/3, 
the shock velocity with respect to the unperturbed gas 
is Vs — 4/3 Vtop, whereas the pressure in the post-shock 
region (i.e., in the cocoon) is expressed in terms of the 
velocity of the piston Vp = Vtop and the upstream gas 
density po as 



3 2 4 2 



(15) 



By comparing this with the pressure balance Ptop — 
Bq/Stt at the contact discontinuity at the top of the 
tower, we see that the tower grows with a velocity of 
order the Alfven speed corresponding to po- 



^top ~ Va 



Bo 



^0 



(16) 



But, as we have shown earlier, Vtop should be of the order 
of _Ro Af2. Thus, we obtain the scaling of Rq with ^fg, Afi, 
and Po: 



(17) 



We can also relate the radius of the tower to the radius of 
the central rotating conductor. Estimating the poloidal 
flux as 'I'o ^ BdR% where Bd and Rd are the typical 
magnetic field and the radius of the the base of the tower 
{e.g., the inner part of the accretion disk), we we get 
-So ^ Bd{Rdl Rq)^ and hence 



Rq 
Rd 



Vd 



1/3 



(18) 



Here Vd = i?t;A57 is the characteristic differential rotation 
velocity of the central conductor and 



Bd 



A,d 



V47rpo 



(19) 



is a composite Alfven speed involving the disk mag- 
netic field and the unperturbed star's plasma density; 
it doesn't have a direct physical meaning and thus can 
be arbitrarily high. 

Notice that, as the tower makes its way through the 
star, ^0 and Af2 remain unchanged, whereas the third 
parameter, the unperturbed density po at the top of the 
tower, changes. It drops rather rapidly for a typical col- 
lapsar progenitor and so the radius of the tower increases 
as the it grows. Thus we expect that, in a realistic sit- 
uation, the tower will not be a straight cylinder, as we 
have assumed here; instead, its radius will be some func- 
tion of the vertical coordinate z. However, as is seen from 
Eq. (fT7| . i?o scales only weakly with po (as Po ^^^), which 
somewhat justifies the constant-radius approximation. 
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The Numbers 

Now let us make some quantitative estimates based 
on the above scaling relationships; unavoidably, these es- 
timates will be very crude. For definiteness, consider 
the accreting black hole scenario for the central engine. 
We assume that the core of the star has collapsed into 
a black hole of fiducial mass M — 3 Mq with a grav- 
itational radius Rg = GMjt? ~ 5 km, and that some 
of the continuously infalling material has formed an ac- 
cretion disk around the black hole. We take a fidu- 
cial disk radius of Rd ~ 6i?g ~ 3 • 10^ cm and an 
initial poloidal flux (per unit toroidal angle) of — 
i?^So ~ lO^^Pisi?^ g 5 in cgs units. This poloidal flux 
is being continuously twisted by the differential rota- 
tion of the disk, with characteristic angular velocity 
Af] = nK{Rd) 4 • lO3sec-i(Af/3M0)i/2 i^-J/^. 

Taking the fiducial stellar background density to be 
po — 10^ g/cm'^, the tower outer radius can be estimated 
as 



-Ro ^ Rd 



VA,d 

Vd 



3Rd 



B 



d.15 



Rd.6.5^^3.5y 



1/3 



(20) 

resulting in Rq ~ 10^ cm for our fiducial parameter val- 
ues. Then we get the following expressions for all the 
other parameters: 



Pn 



Rq \Ro 

niD f Rd.e.5^^3.5\^^^ 1/3 
0.1 Bdl ) Pole 



Bd,15 



10" G: 



(21) 



Va,o 



Bo 



= 3 • 10^*^ cm/sec B^' R^L ^^l'.^ Po.i^^) 



.2/3 



V47rpo 



Notice that our crude estimate results in Vtop ~ Vaa 
being comparable to the speed of light c. Therefore, a 
fuUy-relativistic treatment of the problem would be more 
appropriate (see Sec. for discussion). Such a treatment, 
however, lies beyond the scope of the current paper. 

Also, we can estimate the post-shock pressure in the 
hot cocoon above the tower as 



P, 



top 



B 

-2- ~ 4 • 10^^ erg cm"^ B? 



0.14 ■ 



(23) 



where i?o,i4 = Po/(10^'' G). At such high energy density 
the radiation pressure probably dominates over the gas 
pressure; we can therefore estimate the plasma tempera- 
ture in the post-shock region as 



shock between the cocoon and the unperturbed stel- 
lar material, the baryon density in the cocoon is sim- 
ply Apo, and so the baryon rest-mass energy density is 
4poc^ — 4 • 10^''erg cm~^ po,6: and hence still exceeds 
the radiation/pair energy density by an order of magni- 
tude. The total magnetic energy contained in the tower 
of height Ztop can be estimated as 

B^ 

S,„ag(t) ^ 2nRl Ztop(i) ~ 2-10^" erg R^ j Ztop.g B^ , 

(25) 

which is a noticeable fraction of a typical GRB energy. 



DISCUSSION 

Transition to Relativistic Expansion Regime and the 
Jet Opening Angle 

As we have shown above, the tower radius Rq, and 
hence its growth velocity T4op ^ AflRo scale with back- 
ground density as Pq . Therefore, as the tower ex- 
pands into the outer regions of the star, Vtop inevitably 
reaches the speed of light at some critical density po.rci- 
For our fiducial values Bd = lO^^ G, i?d = 3 • 10^ cm, 
and Afi = 3 • 10"^ sec^^, this critical density for the tran- 
sition to the relativistic regime is po.rci ~ 10® g/cm"^ [see 
eq. ([22]) ]. For a typical massive stellar GRB progenitor 
this corresponds to a distance Z^d from the center on the 
order of 10* cm; according to Eq. (|20|) . the corresponding 
radius of the tower is R^d ^ 10^ cm. 

Our non-relativistic model becomes becomes invalid at 
this point and the subsequent expansion of the tower calls 
for a relativistic generalization of the magnetic tower 
model, which is beyond the scope of this paper. How- 
ever, we probably can derive some physical insight into 
the relativistic regime by looking at the results of fuUy- 
relativistic hydrodynamic simulations '5l'| . In those simu- 
lations the relativistic jet remained coUimatcd as it prop- 
agated through the star, across several orders of magni- 
tude in po- This has been attributed to recoUimation 
shocks in the cocoon and to relativistic beaming in the 
jet. But the physical processes in the cocoon should not 
change if we replace the inner relativistic hydrodynamic 
jet with a relativistic magnetic tower. Furthermore, we 
expect some additional magnetic collimation due to the 
hoop stress. To sum up, we expect the magnetic tower to 
remain coUimated even after it transitions into the rela- 
tivistic regime. In particular, we suggest that the final 
opening angle of the outflow will be about the inverse 
aspect ratio of the tower at relativistic transition: 



top 



top 



v ^ 



1/4 



2 • 10^° K ~ 2MeV, 



(24) 



where a ~ 7.6 • 10 ^ erg cm K . On the other 
hand, since we are dealing with a strong hydrodynamic 



A9 < 



Rq,ic 



0.1. 



(26) 



Whether this prediction is true will have to be deter- 
mined by a fuUy-relativistic analysis and by relativistic 
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MHD simulations, which we hope wiU be completed in 
the near future. 



Small-Scale Structure of the Outflow 

The physical picture presented in this paper, with its 
smooth coherent magnetic structure, is just an idealiza- 
tion used to get the main ideas across in the clearest 
way possible. The actual magnetic field, especially if it 
is produced by a turbulent dynamo in the PNS or in an 
accretion disk, will of course be different from such a sim- 
ple axisymmetric system of nested flux surfaces. Instead, 
it may consist of an ensemble of loops of different sizes 
and orientations. It may thus have a highly-intermittent 
substructure on smaller scales, both temporal and spa- 
tial. However, each of these smaller magnetic structures 
is subject to the same physical processes as the simple 
large-scale configuration: twisting by differential rota- 
tion and a subsequent inflation controlled by the exter- 
nal pressure of the cocoon and of the other loops growing 
at the same time. As a result, a more realistic picture 
may look like a train of spheromak-like plasmoids, push- 
ing each other out along the axis (see Fig. [7]). Hoop 
stress still works inside each of them, and so the overall 
dynamical effect may be qualitatively similar to that of 
a single tower. A similar picture may also develop if a 
large-scale twisted magnetosphere becomes unstable to 
the kink instability and undergoes flux conversion as a 
result, breaking up into smaller plasmoids (see below). 
In either case, the resulting multi-component structure 
of the outflow may be responsible for the observed inter- 
mittency in GRBs. 



Effect of MHD Instabilities on the Tower Evolution 

One of greatest uncertainties in our model is the ef- 
fect of MHD instabilities in the highly twisted magnetic 
structure. In this section we discuss two such instabili- 
ties: Rayleigh-Taylor and kink. 

(i) Rayleigh-Taylor (or its magnetic counter-part, 
Kruskal-Schwarzschild) instability may occur at the in- 
terface between the growing magnetic bubble and the 
overlying colder, denser stellar material. This instability 
is expected to affect purely hydrodynamic fireball models 
as well; if anything, strong magnetic fields are expected 
to suppress it somewhat, although probably not com- 
pletely. The Rayleigh-Taylor instability may cause split- 
ting of a coherent magnetic structure into several sepa- 
rate strands interlaced with stellar matter [sol |40| . For 
example, in the Arons model [iOl, the stellar envelope is 
"shredded" by the nonlinear Rayleigh-Taylor "fingers". 
This leads to creation of several evacuated channels that 
allow the electromagnetic energy produced near the cen- 
tral engine to escape through the star. Arons further 



argued that these channels suffer only a small amount of 
mixing with the non-relativistic stellar material due to 
the Kelvin-Helmholtz instability. In light of his work, we 
think that our magnetic cavity and/or the subsequently- 
formed magnetic towers may also suffer from fragmenta- 
tion into Rayleigh-Taylor "fingers". However, as is well 
known, initially small-scale fingers quickly merge to form 
a few large ones in the nonlinear stage. Therefore, one 
does not expect strong mixing of the baryons from the 
stellar envelope into the magnetosphere. The exact ge- 
ometry of the outflow may be affected somewhat and a 
strong time-variability may develop, but, overall, we ex- 
pect the outflow to survive. More research is needed in 
order to assess the implications of this instability for our 
scenario. 

(ii) As the magnetic configuration is twisted up, it may 
become prone to a non-axisymmetric kink-like instabil- 
ity. This may happen both during the pulsar-in-a-cavity 
phase and during a later magnetic tower phase. Whereas 
the stability of the pulsar-in-a-cavity has not yet been 
studied, several non-relativistic 3D MHD simulations 
have recently addressed [1^ [6l|, H^j the kink instability 
of magnetic towers (although not in the GRB context). 
They seem to indicate that during the first few rotation 
periods, the tower is stabilized by the surrounding high- 
pressure gas, but at later times a large-scale external kink 
does develop. As a result, the tower's general shape be- 
comes helical. This, however, does not immediately lead 
to the overall disruption of the tower; even though the 
configuration is nonaxisymmetric, its main morphologi- 
cal features remain similar to those in the axisymmeitrc 
casejGlj. Similar conclusions have been reached by Naka- 
mura & Meier [HI] in their 3D-MHD study of Poynting- 
flux-dominated jets propagating through a stratified ex- 
ternal medium. In particular, these authors found that 
a steep external pressure gradient forestalls the instabil- 
ity onset. When the instability does eventually develop, 
the resulting helical structures saturate and do not de- 
velop into full MHD turbulence. Important theoretical 
evidence supporting the idea of the external pressure 
stabilization follows from Konigl & Choudhuri's analy- 
sis of a force-free magnetized jet confined by an exter- 
nal pressure [HI]. They showed that a non-axisymmetric 
helical equilibrium state becomes energetically favorable 
(conserving the total magnetic helicity in the jet) only 
when the pressure drops below a certain critical value. If 
this happens and the external kink does go unstable, then 
this non-axisymmetric equilibrium can be interpreted as 
the end point of the non-linear development of the insta- 
bility. 

In addition to the above non-relativistic studies, sev- 
eral first steps have recently been undertaken to under- 
stand the stability of relativistic jets, in particular in the 
framework of relativistic force-free electrodynamics 67 , 
l68| . However, so far as we know, there have been no 
formal stability studies of relativistic magnetic towers 
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to date. Such studies, both analytical and numerical, 
are clearly needed. They would have to take into ac- 
count several stabilizing effects. First, as Tomimatsu 
et al. have found in their linear stability analysis of a 
narrow rotating relativistic force-free jet, rapid field-line 
rotation inhibits the kink instability [6^. Second, we ex- 
pect that the tower expansion should quickly transition 
to the relativistic regime (see Sec. ), eventually reaching 
a very large 7-factor. Once this happens, the relativistic 
time delay may effectively stabilize the tower [l7|. This 
is because MHD instabilities grow on the local Alfven- 
crossing time in the fluid frame and hence much slower 
in the laboratory frame. As a result, even if instabilities 
are excited, they do not have enough time to develop 
before the break-out of the flow from the star. 

If it does develop, the kink is probably the most dan- 
gerous instability and may lead to a significant, although 
perhaps temporary, disruption. Such a disruption, how- 
ever, is not necessarily a bad thing: the tower may be 
able to reform after being disrupted (as is seen in lab- 
oratory experiments 64|) and the resulting non-steady 
evolution may provide a plausible mechanism for rapid 
variability seen in gamma-ray bursts. Also, as a result 
of such disruption, a significant fraction of the toroidal 
magnetic field energy may be dissipated into thermal 
energy [g^ [t^, which may in fact contribute to the ac- 
celeration of the Poynting-flux dominated outflow and 
to powering the prompt gamma-ray emission at later 
times fiol [lil 17, [l3|. In any case, the kink, and espe- 
cially its nonlinear outcome, is a serious issue that needs 
to be addresses in the future. Axisymmetric sausage in- 
stability also needs to be investigated. 

An important aspect of our problem is that we are 
actually interested not so much in the instability onset 
or its early linear development, but rather in the long- 
term (many rotation periods) nonlinear evolution and its 
overall effect on the outflow. An important consideration 
that then needs to be taken into account is the conserva- 
tion of magnetic helicity. Differential rotation leads to a 
continuous injection of helicity into the system (of oppo- 
site signs in the two hemispheres). In its nonlinear stage, 
the kink instability may lead to conversion of some of 
the toroidal magnetic flux to poloidal flux (in our geom- 
etry); however, it will not destroy the magnetic helicity 
accumulated in the cavity. Thus, whatever the resulting 
configuration will be, it will have to be consistent with a 
growing amount of helicity. If some of the new poloidal 
flux becomes detached from the base due to reconnection 
(which may actually be strongly inhibited, see discussion 
in Section ) , then the resulting configuration may resem- 
ble a train of plasmoids (Fig. [7|) , similar to Section . One 
may in fact imagine a cyclic process involving the growth 
of the tower for several rotation periods, followed by flux 
conversion due to the kink, followed by reconnection and 
subsequent ejection of a plasmoid carrying the magnetic 
helicity (and some of the magnetic energy) injected dur- 



ing the given cycle. 

It is also interesting to make the following observation. 
An unbounded relativistic force- free outlow from a rotat- 
ing conductor is expected to be stable. On the other 
hand, in the case of a confined closed magnetosphere 
with field lines subject to differential rotation, such as 
our pulsar-in-a-cavity problem or a magnetic tower, kink 
is expected to develop. At the same time, as we discussed 
in Sec. , the outflow is uncoUimated in the first case but 
is colHmated in the second case. This suggests that there 
may be a deep connection between stability and lack of 
collimation of axisymmetric relativistic force- free flows. 

Finally, we would like to reiterate that a proper treat- 
ment of these problems requires a time-dependent rel- 
ativistic force-free or full (preferably relativistic) MHD 
analysis and simulations (see Sec. ). 



Reconnection 

Another important issue is magnetic reconnection 
across the equatorial current sheet in the pulsar-in-a- 
cavity magnetosphere or across the separatrix current 
sheet in the magnetic tower. This process may, in princi- 
ple, lead to the break-up of a single structure into several 
smaller spheromak-like plasmoids (see Fig. [7]), similar to 
the cyclic bevavior (involving reconnection) suggested for 
magnetospheres of accreting young stars [7l|,|72|- The ex- 
pected size of the plasmoids and their production rate 
are presently not known. 

We would like to remark, however, that fast a recon- 
nection is difficult to achieve deep inside a collapsing 
star[l^. The reason for this is that fast Petschek-like 
reconnection is now believed to be possible only in col- 
lisionless environments, such as the Solar corona. Earth 
magnetosphere, and tokamak plasmas. The plasma in- 
side a collapsing star, on the other hand, is highly 
collisional[15j, so that classical collisional resistivity dom- 
inates over all other non-ideal effects in generalized 
Ohm's law. It is now believed, on the basis of numer- 
ical simulations zS ZJ|, theoretical analysis 75, 7^, and 
laboratory experiments 77| , that reconnection in such a 
situation proceeds in the ver^ ' slow Sweet-ParkerjlilTi 
regime, whereas Petschek's[80j fast mechanism fails. It 
is of course not obvious that this conclusion can be 
extended to the highly relativistic and optically-thick 
electron-positron plasma in the deep interior of a col- 
lapsing star. However, at the very least, this obser- 
vation casts a serious doubt on the possibility of fast 
and efficient large-scale reconnection in this environ- 
ment. As the magnetic tower grows and eventually 
breaks out of the star, however, the plasma cools and 
the particle density in it drops rapidly. At some point, 
the plasma becomes coUisionless (from the reconnec- 
tion point of view) and this opens up the possibility 
of reconnection and the corresponding delayed magnetic 
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energy release. Indeed, post-breakout reconnection in 
relativistic Poynting-flux-dominated outflows has been 
invoked as a plausible mechanism for powering GRB 



emission 



13 0,111, [3, [IS [Bill. 



Nickel Production 

A key issue for long-duration GRBs is the required 
production of ^^Ni. The supernovae observed to accom- 
pany these GRBs (SN-GRBs) belong to Type Ibc 
They are believed to require radioactive ^^Ni to heat the 
ejecta after the initial expansion of the star. The bright- 
est SN-GRBs {e.g., SN1998bw and SN2003dh) require up 
to several 0.1 Mq of ^^Ni, as inferred from peak optical 
brightness, although on average SN-GRBs do_ not need 

56 



more Ni than the local population of SNe[84|. In fact, 

see 



some SN-GRBs {e.g., GRB060505 and GRB060614; 
Refs. [84, 8&,ii|) may produce little or no ^^^Ni [sj. 

In models of core-collapse SNe, ^^Ni is produced in ma- 
terial heated to T > Tn! ~ 5 x 10^ K by the explosion 
shock launched in the core of the star. The amount of 
^^Ni produced depends on the mass inside of the expand- 
ing shock when its temperature declines below Tni. This 
happens when its radius has reached 



( 3E 



V47rar4. 



1/3 



4xl0^£:5icm, (27) 



where E — E^i x 10^^ erg is the explosion energy and a 
is the radiation constant. The mass inside i?Ni depends 
on the progenitor structure and on the expansion or con- 
traction that took place before the shock reached a given 
mass element. In particular, little ^^Ni is produced by 
a shock, even if very powerful, if it is launched into a 
low density material, as it may occur if a weak initial 
explosion pre-inflates the stellar core so that little mass 
remains within a few 10® cm when the subsequent strong 
shock arrives. Production of ~ O.IM© of ^^Ni usually 
requires ~ 10^^ ergs to be deposited isotropically by a 
quasi-spherical shock within ^ 1 sec, so that little pre- 
expansion of the star occurs. Brightest supernovae, e.g., 
SN1998bw, require energies of up to ~ 10^^ ergs to make 
the ~ 0.5Mq inferred from lightcurve modeling. 

The requirement of fast (^ 1 sec) , isotropic deposition 
of energy for production of ^^Ni is a serious challenge for 
models of the SN-GRB central engine, because the GRB 
engine must last for > 10 sec for relativistic ejecta to 
escape the star and because GRBs are believed to be col- 
limated explosions. The high degree of beaming and long 
timescale for energy deposition makes coUapsar jets inca- 
pable ofproducing anywhere near the required amounts 
of ^'^NiQ. Therefore, in the original coUapsar model, 
with a black hole accretion disk as the central engine, 
the ^^Ni is produced by a slower bi-conical disk wind 
that constitutes a distinct explosion component [1, [s^] ■ 



In our magnctar model, on the other hand, ^^Ni can be 
produced behind a roughly spherical shock driven by the 
initial quasi-isotropic expansion of the magnetosphere. 
The expansion becomes coUimated and a magnetic tower 
forms only at a later stage. Thus, the coUimation pro- 
cess involves both a quick isotropic expansion followed by 
a beamed component. We feel that this modification to 
the magnetar scenario, i.e., the inclusion of the magneto- 
sphere interaction with the surrounding star, strengthens 
its viability as a model for the long GRB central engine. 

On the other hand, one cannot exclude the possibility 
that, if most of the energy produced by the central en- 
gine escapes through the narrow jet channel, giving rise 
to a GRB, then there may not be much energy left to 
explode the rest of the star. If this happens, then, even- 
tually, almost all of the stellar material falls into a black 
hole. This would include all the ^^Ni that might have 
been produced, leaving no observable supernova signa- 
ture (e.g., GRB060614). 



Pulsar Kicks 

Since most of the extracted rotational energy of the 
neutron star travels vertically through the two oppositely 
directed channels, a significant amount of linear momen- 
tum is also transported up and down and hence a reactive 
force is exerted on the NS from both sides Even a small 
imbalance in the reactive magnetic force may impact a 
sizable overall momentum to the NS. For example, tak- 
ing the total initial rotational energy of the NS to be 
^-rot = 5 • 10^^ erg, the momentum transported out in 
each direction is P = i?rot/2c ~ lO''^ cgs (it will be even 
larger if the propagation speed is less than c). Thus, just 
a 10% imbalance may result in the NS terminal velocity 
of order of utorm — O.IP/Mns ^ 300 km/sec. 



Suggestions for Numerical Simulations 

In this section we discuss a sequence of numerical stud- 
ies of the interaction between a magnetar's magneto- 
sphere and its birth environment, employing a range of 
plasma descriptions. Each of them will be able to ad- 
dress a subset of key issues with increasing degree of re- 
alism. For example, non-relativistic axisymmetric MHD 
simulations can address the following questions: What 
basic magnetic configuration results when the conditions 
we describe are set up? What is the overall magnetic 
field structure? When do magnetic towers form? How 
does the tower shape change as it expands into lower- 
density regions? How strongly is the magnetic field con- 
centrated towards the axis? How does the Poynting flux 
depend on radius and height? What is the effect of 
the neutrino- or magneto-centrifugally-driven winds? To 
what degree does the cocoon help collimate and stabi- 
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lize the tower? How rapidly do the cocoon waUs spread 
laterally? Can the cocoon expansion result in the disrup- 
tion of the star? In addition, some of the physical pro- 
cesses described in this paper may be relevant to other 
astrophysical systems, including non-relativistic central 
objects {e.g., planetary nebulae 48, H^]). For this rea- 
son, non-relativistic MHD simulations are of interest in 
themselves, as well as a first step toward fully-relativistic 
MHD. For example, recent non-relativistic MHD sim- 
ulations indicate that the magnetic-tower mechanism 
can operate successfully in a variety of astrophysical 
environments H, 57, 58, 0, [8§|, including collapsing 
massive stars [39j. 

Eventually, however, one will have to consider rel- 
ativistic effects outside of the magnetar light cylinder 
or in the inner part of the black hole's accretion disk. 
This can be investigated using the relativistic force-free 
degenerate electrodynamics (FFDE) approach, valid in 
the case of a highly magnetized plasma with negligible 
pressure and inertia. Time-dependent force-free codes 
have recently been successfully used to study pulsar 
magnetospheres[9^ |9l|, |9l|- In the pulsar- in-a-cavity 
context, the cavity wall may be represented by a rigid 
conducting outer boundary. While this case may not be 
directly relevant to the realistic physical environment, 
some basic aspects of a bounded rotating magnetosphere 
may be understood using this description. Furthermore, 
the full magnetar-in-a-star problem can be investigated 
by a hybrid simulation employing a relativistic force-free 
code inside the cavity and a relativistic hydrodynamic 
simulation outside. 

Finally, relativistic-MHD simulations will be able to 
address questions fundamental to the application of 
magnetar-driven magnetic towers to GRBs, including 
the beaming angle and angular distribution of energy 
flux and the growth of instabilities in the relativis- 
tic outflow. There now exist several relativistic MHD 
codes di, [ii, lii, H [93, lii, liil that can be used for 
this problem. Of interest would be a set of simulations 
with a range of plasma (3. The low-/? simulations should 
match on to the FFDE case, at least qualitatively. Once 
these simulations are analyzed and the basic physical pro- 
cesses elucidated, /3 can be gradually increased enabling 
an understanding of how plasma inertia and pressure af- 
fect the dynamics of the magnetosphere expansion and 
collimation. 

The basic scenario described in this paper can initially 
be explored with axisymmetric simulations. However, 
in order to assess the role of non-axisymmctric instabili- 
ties, such as the Rayleigh-Taylor, Kelvin-Helmholtz, and 
kink instabilities (see Sec. ), three-dimensional simula- 
tions will eventually be necessary. One of the goals of 
such an investigation will be to estimate and the degree 
of mixing of the envelope baryonic material into the mag- 
netosphere. In addition, one also would like to study the 
development and interaction of the MRI and the Parker 



instability. 

The operation of the Parker instability, leading to the 
development of a highly-magnetized low-density corona, 
may be strongly influenced by neutrino cooling. Thus, 
a realistic treatment of the neutrino heating and cooling 
processes (as it has been done, e.g., by Burrows et aL[39j) 
is an essential physical ingredient of the overall problem. 

In summary, numerical simulations of the full prob- 
lem, including a detailed description of the central engine 
with relevant microphysical processes and neutrino trans- 
port, are desirable for a comprehensive understanding of 
the formation and evolution of a millisecond-magnetar- 
driven magnetic tower inside a collapsing star. We be- 
lieve that such simulations will very soon become techni- 
cally feasible. 



CONCLUSIONS 

The core collapse of a massive rotating star may result 
in two distinct outcomes, both plausible candidates for 
the GRB central engine. The first one is a stellar-mass 
black hole with an accretion disk. The second is a neu- 
tron star. In this second scenario, the young neutron star 
formed as a result of core collapse has to be a millisecond 
magnetar in order to be relevant for GRBs. 

In this paper we mostly focus on the millisecond- 
magnetar scenario, although many features of our model 
are also applicable to the black-hole case, which we have 
considered in our previous paper (Ref. [la]). Of particu- 
lar interest to us is the interaction between the rapidly- 
rotating magnetar's magnetosphere and the surrounding 
infalling stellar envelope. We argue that the stellar mate- 
rial provides a confining (ram) pressure that has a strong 
effect on both the size and the shape of the magneto- 
sphere. Namely, it can channel the highly-magnetized 
outflow originating from the proto-neutron star into two 
coUimated magnetic towers. 

More specifically, we suggest that the stalled bounce 
shock — a common feature in models of core-collapse 
supernovae — plays a role of a cavity that confines the 
magnetosphere. The cavity's radius, determined by the 
balance between the pressure of the hot neutrino-heated 
gas and the ram pressure of the infalling material, stays 
quasi-stationary at Rq ~ 200 km during the first few hun- 
dreds of milliseconds after the bounce. To get a quali- 
tative physical feeling for what happens to the magne- 
tar magnetosphere during this stage, we introduce an 
idealized fundamental-physics problem that we call the 
Pulsar- in- a- Cavity problem (s^. A large part of our pa- 
per (Sec. ) is devoted to investigating this problem. For 
simplicity, we consider it under the force- free assumption. 
We show that if the radius of the cavity is larger than the 
pulsar light-cylinder radius, the magnetic field inside the 
cavity continuously winds up. Then, the toroidal field 
strength and hence the magnetic spin-down luminosity of 
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the pulsar increase, roughly linearly with time. The mag- from the Keck Fellowship at the Institute for Advanced 
netic energy in the cavity grows quadratically with time. Study. 
We then estimate that in the context of a millisecond 
magnetar inside a collapsing star the magnetic field be- 
comes dynamically important after a few hundred turns. 
This leads to a subsequent revival of the stalled shock and 
may result in a successful magnetically-driven explosion. 
As long as the expansion of the cavity is non-relativistic, 
the toroidal magnetic field inside it remains larger than 
the poloidal magnetic and electric fields. As a result, 
the hoop-stress collimates the Poynting-flux-dominated 
outflow into two oppositely-directed vertical channels. 

We suggest that these magnetic outflows (either in 
the millisecond-magnetar or the accreting black hole sce- 
nario) should be described in terms of the magnetic 
tower concept, introduced by Lynden-Bell in the AGN 
context ^H- Correspondingly, we investigate the prop- 
agation of a magnetic tower inside a star (see also 



Ref. [15|). In our model, we modify Lynden-Bell's pic- 



ture by considering that the tower expansion is super- 
sonic with respect to the unperturbed stellar gas. We 
envision the growing magnetic tower acting as a piston 
that drives a strong shock through the star. The hot 
shocked stellar material between the shock and the tower 
forms a high-pressure cocoon that envelopes the tower 
and provides the coUimating pressure. In other words, 
the tower in our model is confined not by the pressure 
of the background stellar material, but by its inertia; the 
strong shock and the cocoon act as mediators that con- 
vert the inertial support into the pressure support ul- 
timately acting on the tower. The entire configuration 
grows vertically with time and eventually reaches the 
star's surface, thereby providing a narrow baryon-clean 
channel for the Poynting-flux dominated jet, surrounded 
by a less-coUimated hot cocoon outflow. 

Finally, we discuss the astrophysical implications of 
our model for GRBs and core-collapse supernovae, such 
as ^^Ni production and pulsar kicks. In addition, we 
discuss the role of MHD instabilities, most notably, the 
kink, in our scenario. We also assess the prospects for 
magnetic reconnection and find that it should be strongly 
inhibited in the central parts of the collapsing star owing 
to the high plasma coUisionality there. Finally, we outline 
a set of numerical studies that we believe need to be done. 
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FIG. 1: Stalled shock phase of core-collapse explosion. 
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FIG. 2: Schematic time evolution of the main three forces 
responsible for the stalled-shock force balance. 



FIG. 3: Aligned pulsar inside an infinitely-conducting spher- 
ical cavity of radius Ro . The vertical dashed lines represent 
the pulsar's light cylinder of radius J?lc < Ro- After a time 
of order the light- crossing time R) /c, the poloidal field linos 
outside the light cylinder expand somewhat but still remain 
confined within the caA'ity. Because the toroidal magnetic 
field has to vanish at the equatorial midplane due to reflec- 
tion symmetry, the field lines there cannot corotate with the 
star, Os < c/Rb < As a result, differential rotation is 
established in both hemispheres, AQ = f2, — f2s ~ fi* (for 
Rb ^ Rhc), which leads to continuous generation of toroidal 
magnetic flux. 




FIG. 4: At late times, the poloidal magnetic field is pressed 
against the wall by the centrifugal force of the rotating mas- 
sive equatorial sheet. 
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FIG. 5: Development of a magnetic tower in Lynden-Bell's 
(1996) model. Reproduced from Ref. ^l&J by permission of the 
AAS. 




FIG. 6: Main components of a magnetic tower inside. The 
tower grows rapidly and drives a strong shock through the 
star. The shocked stellar gas behind the shock forms a hot 
cocoon whose high pressure confines the tower. Reproduced 
from Ref. fldl l by permission of the AAS. 




FIG. 7: Magnetic tower may have a substructure represented 
by a train of many spheromak-like plasmoids. This situation 
may arise as a result of spatial and temporal intermittency 
at the base of the outflow and/or due to instabilities and 
reconnection in the tower. 



